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Abstract 



We use Omnes representations of the form factors V ,A\ and A2 for exclusive semilep- 
tonic B ^ p decays, and apply them to combine experimental partial branching frac- 
tion information with theoretical calculations of the three form factors to extract | Vub \ ■ 
We find a slightly lower result, \Vub\ = (2.8 ±0.2) x 10^-^ , than the values extracted 
from exclusive semileptonic B ^ n decays, (3.47 ±0.29 ±0.03) x 10"^ [1], (3.36 ± 
Ph; 0.23) X 10^3 [2], (3.38 ±0.35) x lO^^ pj^ ^nd using all other inputs in CKM fits, 

JlLiI (3.55 ±0. 15) X 10^^ [4,5]. The disagreement is greater when we compare to the result 

^ ■ extracted from inclusive 5 ^XJv decays, |y„/,| = (4. 10±0.30exp±0.29th) x 10"^ [5] 

T— I ■ 

> ' 1 Introduction 

in 

^ \ The magnitude of the element Vub of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing ma- 

CN \ trix plays a critical role in testing the consistency of the Standard Model of particle physics and, 

CN ■ in particular, the description of CP violation. Any inconsistency could be a sign of new physics 

^ ■ beyond the standard model. Vub is currently the least well-known element of the CKM matrix and 

O ■ improvement in the precision of its determination is highly desirable and topical. 

• ^ \ \^itb\ can be determined using inclusive or exclusive charmless semileptonic B decays. The in- 
^ ■ elusive method has historically provided a more precise result, but recent experimental [7-12] 
. and theoretical developments [2, 13-22] are allowing the exclusive semileptonic B ^ K method to 
approach the same level of precision. 

Recently [1] we extracted \ Vub \ from combined experimental partial branching fraction information 
and theoretical [lattice QCD (LQCD) and Light cone sum rules (LCSR)] information on exclusive 
semileptonic B it decays. The Omnes representation was employed to provide parametrisations 
of the form factors. The extracted value turned out to be in striking agreement with that extracted 
using all other inputs in CKM fits and in some disagreement with |y„/,| extracted from inclusive 
semileptonic decays. 

The aim of this letter is to extend the above formalism to study the exclusive semileptonic B ^ p 
decay and independently extract \Vub\ from the recent measurements of the partially integrated 
branching fraction by BABAR [8], Belle [9] and CLEO [11, 12]. We will make use of quenched 
LQCD form factor results [23, 24] for the high region, and LCSR values [25] at = 0. Thanks 
to the Omnes representation of the form-factors, we are able to combine all these inputs, as we 
previously showed for 5 ^ ;r decays. 
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2 Fit Procedure 



2.1 Form-factors and differential decay widthi 

The semileptonic decay B'^ p^i^Vi is detemiined by the matrix element of the V — A weak 
current between a B meson and a p meson. The matrix element is 

{p{k,r))\br{l-r5)u\B{p)) = r)*pT^P, (1) 

with form factor decomposition 

lYlg -\- nip 

rug -\- trip q 

where q = p — k is the four-momentum transfer and rj is the p polarisation vector. The meson 
masses are mg = 5279.5 MeV and nip = 775.5 MeV for 5° and p , respectively. In the helic- 
ity basis each of the form factors corresponds to a transition amplitude with definite spin-parity 
quantum numbers in the center of mass frame of the lepton pair. This relates the form factors 
V, Ai and A2 to the total angular momentum and parity quantum numbers of the Bp meson pair, 
= 1^,1^ and 1+, respectively [26]. The physical region for the squared four-momentum trans- 
fer is < < q^six = ('^B ~ '^p)^- If the lepton mass can be ignored (/ = e or /i), the total decay 
rate is given by 

r(5« ^ p-i^y,) = //"^^^^V [Hq')]H\H^q')\' + \H-iq')\' + \H%')\') (3) 

where Gf = 1.16637 x 10^^ GeV^^ is the Fermi constant and X{q^) = {nig + m^ — q^Y — Am^m^p. 
comes from the contribution of the longitudinally polarised p and is given by 

[ml-ml-qA [mB + mp)Ai{q^) - ^'^^]^^ A2{q^) \ (4) 
V '•^ / ^ Mb + nip I 



where k is the momentum of the p in the 5-meson rest frame. correspond to the contribution 
of the transverse polarisations of the vector meson and are given by0 

H^ = -l{ms + mp)M{q^) T ^^^V{q') \ (5) 

The CLEO Collaboration has also measured partial branching fractions of the differential distribu- 
tion [11, 12] 




q^ [X{q^)Y \ 2sm^ OwelH^iq^ 



+ (l-COS0iv£)'|^+(^')P+(l+COS0iv^)2|//-(^2^|2| (6) 

with 9]/^^ the angle between the charged lepton direction in the virtual V7— gauge boson rest frame 
and the virtual W in the 5-meson rest frame. 



^Note a typo in Eq. (1.7) of Ref. [23], the ± sign should be =F, as used in previous papers of the UKQCD Collabo- 
ration [27,28]. 
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2.2 Omnes parametrisations 

We have previously [1,16, 17, 19,20] used a multiply subtracted Omnes dispersion relation [29,30], 
based on unitarity and analyticity properties, to describe B ^ n semileptonic decays. Here, we 
apply these ideas to B ^ p decays and use for (n + l) subtractions [20] 

F{q^) = ^^fl[F{s,){so-sdr'^^"\ a,is)^ H F = VAuA2 (7) 

■^0 ^ i=0 j=OJ^i^' 

where so corresponds to a pole of the form factor F. We fix sq = m\* and sq = 5th = ("^s + ''^p)^ 
for V and Ai and A2 form factors, respectively. In principle, for the axial form factors one should 
use the square of the 1+ 5-meson mass. The mass of this latter hadron is not well established 
yet, but it appears to be heavier than the 15* resonance. Thus and for the purposes of this 
exploratory work, since it would be reasonably far from a/^^^, it is sufficient to employ i'th- The 
parametrisation of Eq. (|7]) amounts to finding an interpolating polynomial for \n[{sQ — q^)F{q^)] 
passing through the points (^•q — Si)F{si). While one could always propose a parametrisation using 
an interpolating polynomial for \n[g{q^)F {q^)] for a suitable function g{q^), the derivation using 
the Omnes representation shows that taking g{q^) = s^ — q^ is physically motivated [20]. 

2.3 Theoretical and experimental inputs 

We have used experimental partial branching fraction data from CLEO [11, 12], Belle [9] and 
BABAR [8]. CLEO and BABAR combine results for neutral and charged 5-meson decays using 
isospin symmetry, while Belle give separate values for B^ — > p~£^Vi and 5+ p^l^Vi decays. 
Belle use three q^ intervals, and we have added in quadrature the two different systematic errors 
quoted for each q^ bin, and combined charged and neutral 5-meson results. We take the resulting 
systematic errors to be fully correlated. BABAR's untagged analysis also uses three q^ bins and 
we have assumed that the quoted percentage systematic errors for the partial branching fractions 
divided by total branching fraction are representative for the partial branching fractions alone and, 
following BABAR, took them to be fully correlated. CLEO determines partial branching fractions 
as a function of both q^ and of cos Q^^l (see Eq. Q) and complete correlation matrices are given 
in [1 1] for both statistical uncertainties and systematic errors that we have used in our fits. 

When computing partial branching fractions, we have used Tgo = l/Lxot = (1.527 ± 0.008) x 
10^^ s [31] for the B^ lifetime. All the branching fraction inputs are listed in Table [T] 

For theoretical form-factor inputs (listed in Table [21), we use the lightcone sumrule (LCSR) results 
at ^2 = of Ref. [25] and lattice QCD resuks from the UKQCD [23] and SPQcdR [24] Collab- 
orations, near q^^^. LQCD inputs have been obtained in the quenched approximation. There is 
therefore an uncontrolled systematic error, which is not fully included in the errors given in Tabled 



2.4 Definition of 

We implement the following fitting procedure. Choose a set of subtraction points spanning the 
physical range to use in the Omnes formula of equation (|7]). Now find the best-fit value of \Vub\ and 
the form factors at the subtraction points to match both theoretical input form factor values and the 
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range 


GeV^ 


cos Qwi range 


10^5^" 




BELLE [9] 


0- 


8 


[-1,1] 


0.62±0.14±0.06 


0.69±0.12 




8- 


16 


[-1,1] 


1.20±0.23±0.11 


1.12±0.15 




> 16 


[-1,1] 


0.53±0.20±0.12 


0.53 ±0.08 


BABAR [8] 


0- 


10 


[-1,1] 


0.73±0.17±0.21 


0.96±0.15 




10- 


15 


[-1,1] 


0.82±0.10±0.13 


0.71 ±0.10 




> 15 


[-1,1] 


0.59±0.07±0.16 


0.68±0.10 


CLEO [11] 


0- 


2 


[-1,1] 


0.45±0.20±0.15 


0.08 ±0.03 




2- 


8 


[-1,1] 


0.96 ±0.20 ±0.29 


0.61 ±0.10 




8- 


16 


[0,1] 


0.75±0.16±0.14 


0.74±0.10 




> 16 


[0,1] 


0.35 ±0.07 ±0.05 


0.39 ±0.06 




> 8 


[-1,0] 


0.42±0.18±0.31 


0.51 ±0.07 



Table 1 Experimental branching fraction inputs for the function defined in Eq. Statistical and 
systematic errors are shown. We also give branching fractions calculated using our fitted form factors and 

\Vub\. 



experimental partial branching fraction inputs. The x function for the fit is 



60 

E 



I I 



C 



-1 
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+ E 

k.l=\ 



5^"-50-"-(|K,|,Fo,Fi,F2) 



5r-5p-"-(|K,|,Fo,Fi,F2) 



^Bkl 



(8) 



where i^;'" are input LCSR or lattice QCD values for y(<5r?),Ai(<5r^) andA2(<?^), and5j," are input ex- 
perimental partial branching fractions. Moreover, F['^"'^^{q^,FQ,F[,F^ stands for each of the form 
factors F = V,Ai,A2 atq^ = qj, and it is given by equation ^ with three subtractions {si,F{si)) at 
(0,Fo), (2(5r^^^/3,Fi) and ('?max,'^2)- The branching fractions 5°™"^*^ are calculated using f'^^"^^^ 
for V,A] and A2 form factors . There are in total 10 fit parameters: V{0), V{2q^^^/3), V(<?max)' 

^2(0), A2(2<5r^ax/3), ^2(^max) ^nd \Vub\- The latter parameter is 



Ai(0), Ai(2^^,j3), Ai(^^J 
used when computing 



We have assumed that the LCSR and LQCD form factor values have independent statistical un- 
certainties and treated the errors listed in Table [2] for the SPQcdR inputs as purely statistical. For 
the UKQCD data we have put the form factor values in the centre of their systematic range and 
use half that range as the systematic error. We have built a covariance matrix where the statistical 
uncertainties (C7,) are uncorrelated and the systematic errors (£,) are fully correlated, leading to a 
60 X 60 covariance matrix with three diagonal blocks. The first 3x3 and second 21x21 blocks 
are for the LCSR and SPQcdR results and have the form Cij = ofdij. The third block is for the 
UKQCD data and has the form Qj = afdij + £/e^ . We will further discuss the effect of the UKQCD 
systematic errors on \Vub \ below. 

The covariance matrix, Cb, for the partial branching fraction inputs is constructed as follows. For 
Belle and BABAR input data, we have assumed independent statistical uncertainties and fully- 
correlated systematic errors leading to an 6 x 6 covariance matrix with two diagonal blocks of 
the form Cg,j = ofdij ± e/£j. For the CLEO input, we use an 5 x 5 covariance matrix C'^J'^ = 

(^if^j'^Bif^^^^^^ + ^i^j'^Bu^^^^^' where we have read off the statistical and systematic correlation 
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[GeV^] 


V 






Ai 




LCSR [25] 





0.323 ±0.029 


0.242 ±0.024 


0.221 ±0.023 


UKQCD [23] 


12.67 


0.684 ±0. 162 


0.439 ±0.067 


+0.000 


0.70 ±0.49 


+0.08 




13.01 


0.714±0.162+j;;§j5 


0.448 ±0.065 


+0.000 


0.71 ±0.46 


+0.08 

\J.\JD 




13.51 


0.763 ±0.1 55 IJJjj 


0.460 ±0.063 


+0.000 

— U.U / J 


0.72 ±0.43 


+0.10 




14.02 


0.818 ±0. 147 IJ^? 


0.472 ±0.059 


+0.000 


0.73 ±0.42 


+0.12 

— U.Ul 




14.52 


0.883 ±0.141 


0.485 ±0.055 


+0.000 

— 070 


0.76 ±0.42 


+0.14 
— 0"^ 




15.03 


0.967 ±0. 137 lj5;5S 


0.498 ±0.051 


+0.000 
n Diss 


0.78 ±0.46 


+0.16 




15.53 


1.057 ±0.1 34 


0.513 ±0.049 


+0.000 

u.uu / 


0.81 ±0.54 


+0.18 
— u.uu 




16.04 


1.164 ±0.15015;^° 


n 529 ± 047 


+0.000 
-0.066 


84 ±0 71 


+0.20 
-0.07 




16.54 


1.296 ±0.1 84 


0.544 ±0.043; 


+0.000 
-0.062 


0.87 ±0.97 


+0.23 
-0.08 




17.05 


1.46 ±0.26 IS 


0.560 ±0.043 


+0.000 
-0.059 


0.90±1.35 


+0.27 
-0.07 




17.55 


1.67 ±0.40 IS 


0.577 ±0.043 


+0.000 
-0.058 


0.90±1.89 


+0.33 
-0.03 




18.17 


2.02 ±0.68 IS 


0.599 ±0.052 


+0.000 
-0.058 


0.9±2.9^ 


-0.4 
-0.1 



SPQcdR [24] 


10.69 


0.51 ±0.26 


0.354±0.085 


0.38 ±0.26 




12.02 


0.61 ±0.28 


0.384±0.087 


0.49 ±0.30 




13.35 


0.74±0.30 


0.421 ±0.089 


0.65 ±0.35 




14.68 


0.93 ±0.31 


0.465 ±0.092 


0.93 ±0.41 




16.01 


1.20±0.32 


0.519±0.097 


1.41 ±0.56 




17.34 


1.61 ±0.33 


0.588±0.108 


2.39 ±1.23 




18.67 


2.26 ±0.55 


0.678 ±0.134 


4.7 ±4.1 



Table 2 Form factor inputs for the function defined in Eq. ([8]l. For UKQCD we show both statistical 
(symmetrized) and systematical enors, while SPQcdR errors include both systematic and statistical uncer- 
tainties (we are indebted with CM. Maynard and F. Mescia for providing us with these form factors). 

matrices stat/sys^ from tables X and XI, respectively, of Ref. [11]. 

We do not consider any correlation between measurements from different experiments, or between 
different sources of theoretical inputs. Nor do we consider correlations between experimental and 
theoretical inputs. 
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3 Results and discussion 



The best fit parameters are 



It/ I 

\yub\ 


= (^Z. /D ± U.zl J X lU 






Vlln^ l'X\ 


— w.Doiinw.w/j 


^(^max) 


= 4.21 ±0.76 


Ai(0) 


= 0.223 ±0.021 


^l(2?max/3) 


= 0.449 ±0.020 


^l(^max) 


= 0.657 ±0.055 


^2(0) 


= 0.231 ±0.022 


^2(2^Lx/3) 


= 0.679 ±0.098 


^2(<?max) 


= 2.76±1.38 



(9) 



The fit has ;f^/d.o.f. = 0.21 for 61 degrees of freedom, while the Gaussian correlation matrix 
can be found in the appendix |Al In figure [T] we show the fitted form factors and the differential 
decay rate calculated from our fit. Partial branching fractions calculated for the same bins as used 
experimentally are given in the last column of Table [TJ Our calculated total branching ratio turns 
out to be (2.30+^26) ^ 10"^, in reasonable agreement with (2.80±0.18±0.16) x 10"^ quoted by 
the Heavy Flavours Averaging Group (HFAG) [31]. 

We have further investigated the effect of the highly asymmetric UKQCD systematic errors on 
\yub\- First, we have completely dropped them and used only the statistical uncertainties on the 
UKQCD points. We find |V„fo| = (2.68 ±0.19) x 10"^. Second, we have performed a Monte Carlo 
where we randomly choose each UKQCD form factor value within its systematic error range, 
with complete correlation between all systematic shifts. For each trial we perform a fit like our 
original one, but setting to zero the systematic errors on the UKQCD inputs. In this case, we find 
\^ub\ — (2.85 ±0.10) X 10^^. Note that this last result is the mean and the standard deviation of 
the fit result for ^ub\ over all the trials, whereas the result above and that quoted in Eq. Q are 
the fit result and error from a single fit. Thus, the result from this second procedure should be 
understood as a shift of +0.09 ± 0. 10 in the value of | V„/,| in Eq. Q. Finally, we have repeated the 
latter procedure, but taking the A2 systematic error to be anticorrelated with those of the V and Ai . 
This results in \Yub\ = (2.86±0.15) x IQ-^. 

From the above discussion, we estimate 

\Vub\=2.S±0.2 (10) 

which constitutes our main result. Quenched approximation systematic effects from LQCD are not 
accounted for by the 0.2 error quoted above. These are difficult to quantify and are a limitation 
here. However, unquenched lattice simulations are now standard and future lattice QCD results 
will address this limitation (although they will also face the problem of an unstable p meson 
for light enough simulated up and down quark masses). Nevertherless, we see that the Omnes 
framework used here provides a fair description of all available experimental and theoretical results 
for semileptonic B ^ p decays, leading to a further independent determination of |Vj,^|. The result 
is lower than the values obtained in the most recent studies of the exclusive semileptonic B 
n decay, (3.47 ± 0.29 ± 0.03) x lO^^ [i], (3.36 ±0.23) x 10"^ [2], (3.38 ±0.35) x lO^^ [3], 
(3.5±0.4th±0.2shape±0.lBR) X 10"^ [21], and using all Other inputs inCKM fits, (3.55 ±0.15) x 
10 [4, 5]. The disagreement is greater when we compare to the most precise result extracted 
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Figure 1 Results obtained from the fit to experimental partial branching fraction data and theoretical form 
factor calculations. The top and the left bottom plots show the three form factors with their 68% CL bands 
(shaded) together with the lattice and LCSR input points (green square LCSR, red dots UKQCD, blue 
triangles SPQcdR). The bottom right plot shows the differential decay rate with 68% CL band (shaded) 
together with the experimental partial branching fractions divided by the appropriate bin-width (histograms 
and points). Green squares, red dots and blue triangles denote BABAR, Belle and CLEO results, respec- 
tively. 

from inclusive 5 ^ XJv decays, \V^b\ = (4.10±0.30exp±0.29th) x 10"^ [gj ^hus, the hints of 
disagreement between inclusive and exclusive/global-CKM-fit determinations are strengthened. 
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A Gaussian correlation matrix 

Here, we give the correlation matrix of fitted parameters corresponding to the best-fit parameters 
in Eq. © 



/ 1 no 


01 

— U.Wl 




n 1 f. 


n in 


n ^1 


9^ 

— U.Zj 


n OS 

W.UJ 


U. 


n 01 \ 




1.00 


0.01 


0.04 


0.00 


0.00 


0.00 


0.00 


0.00 


-0.01 






1.00 


-0.44 


0.05 


0.28 


0.21 


-0.03 


0.16 


0.07 








1.00 


-0.05 


0.09 


-0.08 


0.02 


-0.08 


0.05 










1.00 


0.10 


0.03 


0.17 


0.12 


-0.46 












1.00 


-0.32 


-0.04 


0.34 


-0.30 














1.00 


0.04 


-0.27 


0.47 
















1.00 


-0.04 


0.21 


















1.00 


0.18 


v 


















1.00 / 



(11) 
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